Surface states and their possible role in the superconductivity of MgB2 
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We report layer-Korringa-Kohn-Rostocker (LKKR) calculations for bulk and surface states (SFS) 
as well as the corresponding photoemission intensities of MgB2. Our theor etical results reprod uce 
very well the recent angle resolved photoemission data by Uchiyama et al, sond-mat/011115^ . At 
least two SFS are assigned. Consequences of SFS on the anisotropy of the upper critical fields and 
other properties in the superconducting state of grains in micro-powder samples are discussed. 
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PACS: 74.25.Jb, 74.70, 74.60.Ec, 73.-r 

Recently, the first angle-resolved photoemission spec- 
troscopy (ARPES) measurements on a MgB2 single crys- 
tal were reported by Uchiyama et al. . They observed 
several dispersive states. Among them one has been as- 
cribed to a surface state (SFS). The very existence of SFS 
and subsurface states for MgB2 as well as surface super- 
conductivity have been predicted first in Refs. 2 and 3. 
However, there are significant quantitative differences be- 
tween Ref. 2 and the ARPES data. The effects of SFS 
for the superconductivity in this challenging compound 
have not yet been discussed in detail. In particular, there 
is an unresolved puzzle with respect to the out-of plane 
anisotropy for the upper critical fields = H^^H^2- 
It is moderate for single crystals 1.6 <"fH < 3.5. |j] But 
large values 6 < 7_ff < 9 have been deduced indirectly 
from measurements on pure powder samples with small 
grains. Enhanced superconductivity due to locally 

enhanced DOS (density of states) has been predicted ||] 
for idealized B- as well as Mg-terminated surfaces (BTS, 
MgTS) although in conflict with experimental data. We 
report on LKKR calculations for MgTS and BTS and 
compare them with the ARPES data. ||^ Furthermore 
we propose a scenario which involves essentially SFS to 
resolve some puzzles of general interest mentioned above. 

The electronic DOS for the semi-infinite system MgB2 
as well as the photoemission spectra, were calculated us- 
ing a LKKR code 1^ . This code treats the photoemission 
process within the "one-step" approach. Dealing with 
the Green function formalism, it is capable to work with 
a complex potential, allowing the description of broad- 
ening effects due to the finite lifetime of quasi-particles. 
However, a set of input parameters must be provided 
phenomenologically in order to match experimental data 
as close as possible. As an input to the LKKR code, 
the real part of the one-electron effective complex po- 
tential was calculated self-consistently with the help of a 
LMTO (linear muffin-tin orbital) code. This preliminary 
calculation was done for the infinite MgB2 crystal with 
a basal plane lattice constant of a=3.085 A and a ratio 
c/a=1.142. The von Barth-Hedin correlation energy 
functional was used in connection with the local den- 
sity approximation (LDA). The resulting Fermi energy 
Ep « 13.1 eV measured from the muffin-tin zero, while 



the muffin-tin radii were 1.711 A and 0.988 A for the Mg 
and B, respectively. A direct evaluation of the self-energy 
correction to be added to the LMTO potential is noto- 
riously very difficult. Hence, we fixed it empirically. As 
for its imaginary part we chose different values for states 
below and above the Fermi level, -0.05 eV and -1 eV 
respectively. Its real part was chosen to be -1 eV and ap- 
plied to the photo-electron final states only. Rigorously, 
the full complex self-energy correction is not constant in 
the energy intervals under consideration, but here such 
details can be omitted without losing much accuracy. 

The geometry of the semi-infinite MgB2 system was 
simply modeled cutting the ideal crystal perpendicularly 
to the c axis. Hence, no possible relaxation or recon- 
struction effects were taken into account. The topmost 
surface layer consists either of Mg atoms or B atoms. In 
the MgTS-case a work function cj) —4.2 eV was used, in 
the BTS =6.1 eV, as suggested in Ref. §. In both 
cases the surface potential barrier was modeled as a step 
potential, reflecting for occupied states and non reflect- 
ing, but refracting, for unoccupied ones. Its position d 
was fixed to 0.987 A with respect to the topmost layer, 
while its height is given by the sum Ep + 4>. The energy 
of impinging photons at 28 eV and the linear polariza- 
tion _L to the symmetry directions S and A were chosen 
in accordance with the experiment. |0| 

In Fig. |l| we show the DOS of the surface layer in 

DOS contribution ot Mg in ttie surtace layer 
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FIG. 1. The DOS in the topmost layer (MgTS-case) vs 
binding energy and ky-vector along the fM (E) and f K (A) 
symmetry directions. The DOS is measured by the darkness 
on a logarithmic grey-scale. The step-like surface potential 
barrier is at a distance d=0.99 A from the topmost layer. 
Notice the Mg-derived SFS at eV at f . 
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the case of a MgTS. The quasi parabohc feature with 
a minimum around -0.5 eV on the T point is clearly a 
SFS, lying in a gap in the fc^-projected bulk band struc- 
ture 0. It is s, Pz, in character, in a proportion 
roughly 3:2:1 at the f point. Its characteristic decay 
length is 2.6 A, whereas the distance between the Mg 
and B planes is 1.762 A. In Fig. ^ the dispersion of bulk 
bands is also visible. The narrower stripes correspond 
to the quasi 2-dimensional a bonded Pa;,y-bands, while 
the broad parabolic tape corresponds to the tt bonded 
p^-band. 

SFS energies are much more sensitive to the values of 
the surface potential barrier positions d than to the work 
functions (f>. Placing the barrier at a larger d value results 
in a downward energy shift (cf. Figs. I and |). Thus, the 
situation depicted in Ref. 2, where the energy of the SFS 
of a MgTS is « -2 eV at F, can be achieved with a surface 
potential barrier placed at c?=1.32 A above the topmost 
Mg layer. A different situation occurs in the BTS-case 
(cf. Fig. ^). Here the result of Ref. 2 is reproduced by a 
smaller d-value than our empirical one. This is probably 
due to the in-plane B-B distance which is 57% smaller 
than the Mg-Mg one, so that in a BTS slab calculation, 
the image plane position could have been chosen nearer to 
the surface than in the MgTS. The SFS which we found 
for a BTS in the considered energy range, runs parallel 
to the Pz broad band near its upper edge. It is depicted 
by the solid line peak at a binding energy « -0.8 eV in 
Fig. H calculated at ky = 0.7-M. It decays towards the 
bulk with a characteristic length of 9.8 A and exhibits 
Pz symmetry. If steps are present on the real surface, 
yielding islands of Mg- and BTS, a single effective 
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FIG. 2. The DOS of the surface layer in the case of a 
MgTS as a function of the binding energy, for ky — T. The 
different curves were calculated with step-like surface poten- 
tial barrier placed at various distances d (given in A) from 
the surface. The imaginary part of the self-energy was set to 
Iml!=-10 meV in order to sharpen the peaks. 
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FIG. 3. The DOS of the surface layer in the BTS case 
as a function of the binding energy for ky = 0.7-M. The 
curves were calculated using a step-like surface potential bar- 
rier placed at different distances d (given in A) fr om the top- 
most layer. ImS as in Fig. 2. 

(averaged) d value could be justified to analyze the 
ARPES data. In view to have either a MgTS a BTS (or 
even both if we have a step- like surface), we calculated 
the ARPES spectra shown in Figs. 4 and 6. 

We would like to stress that the crucial effective d pa- 
rameter must be fitted from experimental data, and the 
only available data jl] is shown in Fig. ^. Nice agree- 
ment with our calculations is obtained, if one superposes 
Figs. 4 and 6, suggesting that the measured sample pre- 
sented both Mg- and BTS. According to our calculations 
the feature centered at «0.4 eV at f is then to be as- 
signed to a TT-bonded Mg derived SFS in the case of 
MgTS. The parabolic feature centered at «2.7 eV is most 
likely, in our opinion, to be assigned to a 7r-bonded B 2pz- 
derived SFS arising from a BTS, in accordance with the 
calculations in Ref. 2. However, we admit that from the 
presented experimental data, an alternative assignment 
of this parabolic feature to the bulk B 2p2 derived 7r-band 
edge cannot be completely excluded. An analogous con- 
sideration holds for the B a derived narrow bulk stripes, 
where a surface resonance state is found with our d pa- 
rameter value (cf. Fig. 3). Further ARPES experiments 
at variable photon energies should be performed to solve 
the full assignment problem. We emphasize that our in- 
terpretation mainly in terms of SFS does not reduce the 
importance of Ref. 1 which is the first experimental veri- 
fication of the correct description of the electronic struc- 
ture of MgB2 by the LDA. This is very important with 
respect to some times discussed strong many-body ef- 
fects (correlation or non-adiabaticity) which would cause 
pO| , pl| ) strong renormalizations of bulk and SFS as well. 
If strong correlation effects were present, they should 
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FIG. 4. Theoretical ARPES energy distribution curves 
taken along E (left panel) and A (right panel) symmetry direc- 
tions for a MgTS. The intensity is measured by the darkness. 
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FIG. 5. ((Color), after Uchiyama et al.[l\) Experimental 
ARPES data (upper panel) . The intensity is measured by the 
brightness obtained from the second derivative of the ARPES 
spectra. The assignment of one SFS (red) and bulk states (vr 
(blue) and a (green)) proposed in Ref. 1 (lower panel). 
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FIG. 6. The same as in Fig. 4 for a BTS. 

be observed first in the SFS, because the screening of 
Coulomb interactions is weaker near the surface rather 
than in the bulk. Taking into account the perfect agree- 
ment of LDA with the ARPES data, we conclude that 
such many-body effects should be of minor importance. 
With respect to surface superconductivity to be discussed 
below, the experimental determination of the electron- 
phonon (el-ph) coupling constant A from the temperature 



dependence of the ARPES spectra like in Be would be of 
great interest |l^ , |l3|] . 

Now we address possible consequences of the SFS with 
respect to physical properties of various samples with em- 
phasis on superconductivity and micro-powder samples. 
Since for typical samples the size in c-direction is about 
10 nm only, about one quarter of the volume is affected 
by the presence of SFS. Since for MgTS the SFS shown 
in Fig. 1 is derived from orbitals with cylindrical symme- 
try and there is no contribution from orbitals with px^y 
character, a very weak local el-ph interaction is expected. 
Thus it acts as a bad or non- superconducting layer in 
the sense of the proximity effect. Due to the interac- 
tion with the bulk states the resulting superconductivity 
near a MgTS will be weakened even for a clean MgTS, 
at variance with the view put forward in Ref. 3. 

Thus the presence of such SFS might explain the ab- 
sence of a third upper critical field i?c3 and the local 
suppression of superconductivity in relatively weak exter- 
nal magnetic field slightly above 2 Tesla This was 
regarded as evidence for a small value of i?^2i i-e. a big 
anisotropy 7// ~ 6 to 9. Anyhow, this would be in sharp 
contrast to the anisotropy caused by the Fermi velocities 
7^ ~ 6 to 7 which yields 7|^ ~4.2 to 5.5 as an upper 
limit iQ derived for the strongly interacting subsystem 
of (T-holes. [|l^ Due to the two-band character of MgB2 
]l5|-p7| which involves also the weakly interacting nearly 
isotropic 7r-electrons the total 7// is further reduced till 
below 4 to 5, provided any order parameter anisotropy of 
the two different gaps can be neglected |l^,|l^. A mod- 
erate jH 3 is in accord with available single crystal 
data. Q 

The interplay of superconducting bulk and SFS sub- 
systems considered here contrasts the approach adopted 
in Ref. p9) where due to symmetry breaking at the sur- 
face mixing and pair breaking of a- and tt- states in the 
sense of a two-band model occurs. Notably, the absence 
or presence of superconducting SFS may affect the sur- 
face pinning. Bulk vortices will be attracted or repelled, 
respectively. This might explain the strong variations of 
pinning properties after applying pressure, which might 
affect the interfaces of the grains. Furthermore the pres- 
ence of surface superconductivity on SFS will affect the 
tunneling data. In the context of frequently considered 
3D multi-band models for MgB2 it is important to make 
a correct assignment of the observed superconducting 
gaps. In particular, in view of the ARPES interpretation 
we suggested above, the gap at about 5 meV observed 
by high-resolution (angle integrated) photoemission spec- 
troscopy po| might be ascribed to superconductivity on 
SFS. 

Up to now two-dimensional (2D) electrons in metals 
were observed only in artificial layered structures plf . 
Let us consider in short the new physics which can be 
gained by studying generic 2D electrons at the surface of 
a 3D superconductor. As most interesting we consider 
the possibility for enhancement of the superconductivity 
of these 2D electrons due to increased local DOS or soft- 
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ening of local phonon modes at the surface. In this case 
in the framework of Ginzburg-Landau theory, for exam- 
ple, we have to introduce a 2D superconducting order pa- 
rameter ^21) in parallel to the bulk one ^'au; additionally 
the phenomenology could be complicated by necessity of 
taking into account pairing in different and weakly in- 
teracting bands. As a whole the surface enhancement of 
superconductivity could remind the twinning plane su- 
perconductivity observed in Sn and Nb. In this case 
the dependence of magnetization as a function of temper- 
ature and magnetic field can be similar to the fluctuation 
conductivity of the bulk. Roughly speaking, due to the 
proximity effect 2D conductivity can create a "frozen" 3D 
fluctuation. A systematic study of fluctuation magnetism 
in powdered samples with different grain sizes can give 
important hints in this direction. It will be very infor- 
mative to perform investigation of fluctuation magnetism 
psf and to determine the mass anisotropy of fluctuation 
Cooper pairs using MgB2 microcrystals sintered under 
high pressure p^ . 

However it is more plausible that the influence of the 
surface is to reduce the local superconductivity: we have 
to take into account despite the possible Mg or B ter- 
minations of the crystal, also the adsorption of various 
impurities etc In this case the influence of the sur- 
face can be traced out in the vortex pinning in the super- 
conducting state. Indeed, strong asymmetric hysteresis 
loops were observed after the grinding of a bulk sample 
into fine powder |25|. Surface pinning was reported for 
a powder sample pressed into a pellet without sintering. 
p6[ Both possibilities, enhancement and suppression of 
superconductivity yield additional surface pinning. Fi- 
nally, we argue that 2D electronic states can give an im- 
portant contribution to the current dependent correction 
of the work function of the superconducting metal, the 
so called Bernoulli potential. 

We have shown that the available ARPES data can be 
well described by the LDA. In particular, there is clear 
evidence for the presence of at least two SFS. Since the 
SFS depend very sensitively on the nature of the ter- 
minated layer, comparing our calculated intensities with 
the measured ones, we suggest that the studied single 
crystal contained both types of terminated layers. There- 
fore the investigation of especially prepared surfaces with 
only one predominant surface layer is of interest. In this 
context the theoretical study of steps between B and Mg 
terminated layers would be interesting. Since the SFS are 
generated by the large gapped region in the BZ, similar 
phenomena might be expected for other layered super- 
conductors such as NbSe2 and even heavy fermion com- 
pounds URu2Si2 or borocarbides where also no Hcs 
has been observed. Taking into account all those possibil- 
ities mentioned above one can expect that comprehensive 
studies of MgB2 could generate interesting physics of the 
surface/ interface of superconductors. 
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